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Abstract 



(N 
> 

CN , We study numerically and analytically the spectrum of incidence 

matrices of random labeled graphs on N vertices : any pair of vertices 
("^ ■ is connected by an edge with probability p. We give two algorithms to 

^P . compute the moments of the eigenvalue distribution as explicit poly- 

nomials in N and p. For large A^ and fixed p the spectrum contains 
C^ ■ a large eigenvalue at Np and a semicircle of "small" eigenvalues. For 

large N and fixed average connectivity pN (dilute or sparse random 

matrices limit) we show that the spectrum always contains a discrete 

^ ' component. An anomaly in the spectrum near eigenvalue for con- 

O ■ nectivity close to e is observed. We develop recursion relations to 

compute the moments as explicit polynomials in pN. Their growth is 
slow enough so that they determine the spectrum. The extension of 
^ ■ our methods to the Laplacian matrix is given in Appendix. 

Cd I Keywords: random graphs, random matrices, sparse matrices, inci- 

dence matrices spectrum, moments 

1 Introduction 

The spectral properties of the incidence matrix of random graphs have mo- 
tivated a large number of studies over the last decades. The same problem 
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is described under rather different names, depending upon tlie aspects tliat 
are under focus and tlie metliod of attack. 

Tlie interest in tliis problem lias several roots in physics. The replace- 
ment of complicated hamiltonians by large random matrices has proved very 
efficient in the analysis of the spectral properties (culminating in level spac- 
ing distributions) of large nuclei since the pioneering works of Wigner and 
Dyson. For many properties, the details of the probabilistic laws govern- 
ing the distribution of matrix elements are irrelevant, and there is a very 
powerful notion of universality. Further motivation for considering precisely 
random incidence matrices comes from several systems in condensed matter 
physics, a good example being conductors with impurities. The pure system 
is modeled by a lattice, and electrons can move along bonds. Impurities 
break bonds. So the hamiltonian can be approximated by the incidence ma- 
trix of the lattice with random bonds removed. If one considers several large 
samples differing by the impurity concentration c, the following properties 
are observed. When c is large, only small islands of metallic atoms exist. 
If c decreases to reach a certain threshold, a large island of metallic atoms 
invades the system. This is classical percolation. However it is generally 
believed that the system remains insulating (the wave function of the elec- 
trons are all localized) until another threshold in the impurity concentration. 
Then some delocalized states appear and the sample is a conductor. This 
is called quantum percolation, a kind of Anderson (de) localization. The an- 
alytic study of this problem on a 3d lattice with random bonds removed is 
very difficult, and this motivated people to look at the much simpler problem 
of a random graph. This forgets about the spatial structure and is a kind of 
mean field approximation. 

In the random graph model, the lattice is replaced by the complete graph 
on A^ points: any two points are connected by a bond (by an edge in the 
language of graph theory). Then, bonds are randomly removed, leading to 
a random graph where only a fraction p of the initial bonds remains. In 
the simplest case, bonds are removed with probability 1—p independently of 
each other. The model can be made more complicated by choosing randomly 
a sign for each bond present in the random graph. This allows interferences 
if the probability amplitude for an electron moving on the random graph is 
the product of the signs of the visited bonds. 

The topology of large random graphs was investigated about four decades 
ago by Erdos and Renyi ||Erdos6(]|] in a remarkable series of papers. The idea 



is to let p vary with A^. There are quite a few different regimes. The most 
relevant for further physical investigations are: 

• The edge-probability p remains fixed as A^ goes to infinity. 
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• The average connectivity a = pN remains fixed as A^ goes to infinity. 

In the first case, the infinite random graph is connected and in a precise 
sense two infinite random graphs of given p are isomorphic with probabihty 
1. The second case exhibits a percolation transition at a = 1. For small a all 
connected components are finite, and only trees contribute to the extensive 
quantities. But for a > 1, a finite fraction of the points lies in a single 
connected component. 

As explained above, the spectral properties of the (signed) incidence ma- 
trix of the random graph, a symmetric matrix with 0, (±)1 matrix elements 
have a great physical interest. Quite often, authors concentrate on the case 
when the distribution of random signs is symmetric. Numerical simulations 
and analytic (mostly supersymmetric) methods have given a great amount 
of information. 

For large A^, fixed p and symmetric random signs, it is known that the 
spectral distribution is a semicircle ||Rodgers"88l and that the level correla- 



tions are those of the Gaussian orthogonal ensemble (GOE), one of the four 
standard universality classes governing random spectra ||Mirlin9l ]. 



The finite connectivity limit has also attracted a lot of attention, under 
the names of dilute or sparse random matrices. It has been argued |[Mirlin9l 



Evangelou92| , |Evangelou92a|| that there is some value a^ > 1 for which de- 
localized eigenstates appear. So this simple model is believed to exhibit a 
quantum percolation transition. 

Our aim in this paper is to use combinatorial methods to compute ex- 
plicitly moments of the spectral distribution for given A^ and p, and in the 
finite connectivity limit. 

Section 2 gives a formal definition of the model and recalls some of its 
topological properties. Section 3 gives the enumerative algorithm for mo- 
ments as polynomials in A^ and p. Section 4 concentrates on the fixed p large 
A^ limit, first numerically (spectrum and level spacing) and then analytically. 
Section 5 deals with the finite connectivity large A^ limit, starting with nu- 
merical computations. In particular, we observe a quantitative change in the 
spectrum near the eigenvalue A = for a ~ 2.7. We give qualitative argu- 
ments for the presence, location and size of delta peaks in the spectrum. Then 
we derive a formal expression for the moment generating function and give a 
recursion relation for the moments, that we use to control their growth. In an 
Appendix, we show how our algorithms can be extended when the incidence 
matrix is replaced by the Laplacian matrix of a random graph. 
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2 The model 

For A^ = 1, 2, ■ ■ -, we define the sample space Q^ as follows : the elements of 
Qn are the symmetric N x N matrices M = (Mjj)j ,,g[i^^], with Mj^j = on 
the diagonal and Mj .,• = or 1 for z 7^ j. So Q^ is a discrete space consisting 
of 2^(^~^)/^ points. We observe that fi^r is in one to one correspondence with 
the set of labeled simple graphs on A^ vertices : to M G Qn, "we associate 
the graph with vertex set {1, ■ ■ ■ A^} and edge set E{M) = {{i,j} \ Mij = 
1}. Then \E{M)\ = l^ij Mij = Y.i<j M^j is just the number of edges of 
the graph associated to M. The word simple above refers to the fact that 
the graphs we consider do not contain multiple edges or edges with only 
one vertex. In the sequel, graph always means simple graph, and we talk 
indiscriminately of matrices or associated graphs. 

For any p G [0, 1] we turn fi^r into a probability space : the weight of M G 
fijv is P(M) = p^(^^)(l -p)mN-i)/2-E(M)^ rp^ rephrase this formal definition, 
the entries of M above the main diagonal are independent random variables 
with the same Bernoulli (binomial) distribution: for i < j, Mij = 1 (or 
equivalently the vertices i and j are connected by an edge) with probability 
p and Mij = with probability 1 — p- 

The quantity a = pN appears as the average connectivity]^, i.e. the av- 
erage number of vertices j connected by an edge {i, j} to a given vertex i. Re- 
mark that this connectivity fluctuates, in contrast to regular graphs ||Jakobson99 | 
where the connectivity is fixed for all vertices. 

We can define a variant of this model by introducing a random sign, 
with a parameter a G [0,1]: for i < j, Mij = +1 with probability ap, 
Mij = — 1 with probability (1 — a)p and Mij = with probability 1 — p. 
The even model a = 1/2, which gives (Mij) = 0, has been studied by some 
authors [[Rodgers88| , [Rodger s90| , [Mirlin91| , |Evangelou92| , Pilvangelou92a)| . If a 



random graph contains no loops (i.e. closed circuits), the parameter a is 
not relevant because all the negative signs can be changed in positive ones 
by a simple change of basis. More generally, it is true if the graph has no 
"frustrated" loops, i.e. no loops with odd number of negative edges. We 
will see later that the random spectrum is not sensitive to a in the large N 
limit with fixed a. So, without explicit indications, we will speak about the 
signless model defined previously, which has a = 1. 

If X is any random variable on Qj^, we use the notation X for the expec- 
tation value (or average) of X. In the next sections, we shall be interested 
in the asymptotic behavior of the moments of the spectral density Tr M ^ 

"'^The terminology average connectivity seems to be well established in the physics liter- 
ature, and we stick to it. It would however be more consistent with general graph theory 
to call it the average degree (of vertices) 
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when A^ ^ oo, first with p fixed and then with a = Np fixed. But first, we 
recall a few fundamental facts on the topology of random graphs. The basic 
reference is [[Erdos60|| ; for a textbook presentation and more references, see 



e.g. Pollobas98 



It is well established that this model has a percolation transition at a = 
1. In the regime a < 1, with probability 1 in the large N limit, all the 
connected components of the random graph are finite: moreover they are 
mostly trees, there is only a finite number of one loop components and no 
other connected components. For a given tree T on n vertices, the average 
number of connected components isomorphic to T is 

^ a"-^e-'^" + o{N) (1) 



\Aut{T)\ 

where |/lMt(G)| for a given graph G is defined as the order of its automor- 
phism group, formed by the permutations of the vertices that leave its inci- 
dence matrix invariant. 

In the regime a > 1, Eq (^ remains valid, but one "giant" connected 
component (equivalent of the percolation cluster for regular lattices) appears, 
with a finite fraction of the N vertices and many loops. This fraction is an 
increasing function of a, covering [0, 1] when a runs from 1 to cx). In the limit 
N large with p fixed, the random graphs is made only with one component. 

The percolation transition at a = 1 is of second order. As usual, critical 
exponents can be defined: by example, the biggest component has a size of 
order iV^/^. 

Moreover this model exhibits ||Mirlin9i| an Anderson localization transi- 



tion, also called quantum percolation transition, at a value a^ > 1. In the 
phase a < aq, all eigenvectors of the random incidence matrix are localized. 
On the other hand, for a > aq, the eigenvectors for which the absolute value 
of the energy is below a threshold E{a) becomes extended. 

By studying nearest level spacing between eigenvalues, we expect an ex- 
ponential distribution of spacing in the localized phase (because the spa- 
tial covering between different eigenvectors vanishes), and a GOE distribu- 
tion ||Mehta9l| in the delocalized phase. With this kind of criterion, the loca- 



tion transition has been numerically estimated [[Evangelou92| , [Evangelou92 



at aq ~ 1.4. By considering quantum percolation on a randomly diluted 
Cayley (or Bethe) tree ||Evangelou92| , [Evangelou924 [Harris82|] , it has been 



conjectured that aq is given by aq loga^ = 1/2, (leading to aq ~ 1.4215299). 
We argued ||BauerOO"a[] that this value is exact for random incidence matrices. 



but that loops have nevertheless some influence on the localization properties. 

So the percolation transition has a drastic effect on the topology of the 

random graph and the localization transition changes the behavior of eigen- 
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vectors of its incidence matrix, but as we shall see later, the transitions have 
no obvious impact on the moments of the spectral distribution of a random 
matrix. In fact, the situation may look paradoxical: many relevant quantities 
for the spectrum of random incidence matrices (for instance the moments) 
can be computed by looking at local structures on the random graph. For 
fixed a and large A^, such structures are trees with probability 1. Hence loops 
that appear for a > 1 seem to play no role. However, the presence of loops is 
crucial to ensure that the statistics of finite structures varies smoothly with 
a. For example, if, instead of random graphs, one considers random forests 
(union of trees) |[BauerU0"5| , one finds that for a < 1, the thermodynamical 



properties are exactly equal to the ones for the random graph model, but the 
transition at a = 1 (when an infinite tree appears) changes the distribution 
of local structures, and for instance the moments of the spectral distribution 
are not analytic at a = 1 for random forests. 

3 Computation of moments 

In this section, we derive Eq. (0), valid for any A^ and edge parameter p, 
which allows, for a given k, to compute directly Tr M'^, (i.e. A^ times the 
kth moment of the density of states), when M is a random incidence matrix 
in Qjsf. A sum rule for p = 1 is given. Then the algorithm is adapted for the 
variant of the model with random signs. Finally we give a compact formula, 
Eq. (^, for the generating function of moments. 

3.1 Direct computation 

For any random incidence matrix M in Q^, Tr M^ = N and Tr M^ = 0, so we 
may assume that k > 2. By definition, Tr M'^ = Y,i[ ... i,.=i Mi^i^Mi^i^ ■ ■ ■ Mi^i^ . 
Because the diagonal matrix elements of M vanish, we can restrict the 
above sum and consider only fc-plet (ii, ^2, ■ " " ? ^fc) such that ii ^ 22, "^2 7^ 
23, ■ ■ ■ , ik-i 7^ ik, iki^ii- We call such /c-plets admissible. 

So, start with an admissible A;-plet / = {11,12^ ' ' i^k)- To compute 
Mi^i^Mi,^i.^ ■ ■ ■ Mi^i^ we argue as follows : the product Mi^i^Mi^i^ ■ ■ ■ Mi^i^ can 
take only two values, or 1. It has value 1 if and only if each factor has value 
1, that is if and only if {ii, 22}, {^2, "^3}, ■ ■ ■ , {4, "^i} are edges of the graph with 
incidence matrix M. From our definition of probabilities on the space of in- 
cidence matrices, this happens with probability p\ where / is the number of 
distinct pairs among {zi, 12}^ {12, h}i " " " ) {^i ^i}- Obviously, I depends on I. 

If we could find an efficient way to count the number of admissible fc-plets 
/ with given /, the problem would be solved. We have not been able to do 
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so. However, there is a simple way, which we now expose, to group together 
famihes of admissible /c-plets that are guaranteed to have the same /. 

Fix a /c-plet W = (f i, ■ ■ ■ , Vk) of elements of {1, ■ ■ ■ , k} with the following 
properties : i) vi j^ V2, ^2 7^ f 3, ■ ■ ■ , f fc-i 7^ Vk, Vk 7^ Vi, and ii) if Vjs > 1, there 
is a. P' < f3 such that vp/ = Vf^ — 1. Such a A;-plet will be called a normalized 
/c-plet in the sequel. 

The first condition is almost the definition of an admissible /c-plet, the 
only difference being that the members are in {1, ■ ■ • , A;}, not {1, ■ ■ ■ , A^}. The 
second condition says that the order of appearance of elements of {1, ■ ■ ■ , k} 
in the sequence (vi, • ■ ■ ,Vk) is the natural order. Because of these two con- 
ditions, wi = 1 and f 2 = 2 in any normalized /c-plet, but ^3 could be 1 or 
3. 

By condition ii), the integers appearing in W build a set of the form 
V = {!,■ ■ ■ ,n} for a certain n < k. Let E be the set whose elements are the 
(distinct) pairs among {vi, V2}, {v2, V3}, ■ ■ ■ , {vk, vi}. 

Now choose an injective map a from V to {1, ■ ■ ■ , A^}. The number of 
such maps is N- = N{N — 1) ■ ■ ■ (A^ — n + 1). Then the sequence / = 
(o"(wi), ■ ■ ■ , <7{vk)) is an admissible /c-plet by the injectivity of a and property 
i) of the sequence W. Moreover, for the same reasons, the number / of distinct 
pairs among the k pairs {a{vi),a{v2)},{a{v2),(T{vs)}r ■ ■ ,Wi'^k),<^{vi)} is 
exactly \E\, the number of elements of E. This number depends on W, but 
not on a. 

More precisely, there is a one to one correspondence between admissible 
/c-plets / and pairs (W, a). Note that the source of a depends on W. 

The bijection involves a simple but useful general algorithm, which we 
call the "label and substitute algorithm". We use it several times in the 
sequel. If (O, O', O", ■ ■ ■) is any finite or infinite list of items (some items can 
be repeated), one can label the items in order of first appearance. This means 
that the first item receives label 1, then the next item different from the first 
one receives label 2 and so on. This gives a one to one map, the "labeling 
map". After that, by replacing each item in the list by its label, one obtains 
a sequence of integers, the "substitution sequence" . It has the property that 
the first term is 1 and that if integer i + 1 > 2 appears as a term, then integer 
i has appeared before. Note that this new sequence is invariant if we apply 
the "label and substitute algorithm" to it. Note also that the knowledge of 
the "labeling map" and the "substitution sequence" allows to reconstruct the 
original sequence. Let us give an example. The list (eat, work, eat, sleep, 
work, eat, work, sleep) leads to the "labeling map" (eat -^ l,work -^ 2, sleep 
— > 3), and to the "substitution sequence" (1, 2, 1, 3, 2, 1, 2, 3). 

If / = {ii,i2,- ' ' y'^'k) is an admissible /c-plet, we apply the "label and 
substitute algorithm" . Then a is the inverse of the "labeling map" of / and 
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W is the "substitution sequence" of /. The properties of I making it an 
admissible A;-plet and the definition of the "label and substitute algorithm" 
ensure that 14^ is a normalized /c-plet. 
Written symbolically, this means that 

E = EE 

I W cr 

where the sum over I is over admissible /c-plets, the sum over W is over 
normalized A;-plets, and the sum over a is over injective maps as described 
above. Inserting Mi^i^Mi^i^ ■ ■ ■ Mi^i^ on both sides of this identity, we obtain 
our first important formula: 

TWi^=^ArMpl^l, (2) 

w 

where V and E are functions of W as defined above. In this formula, the A^ 
and ^-dependence are completely explicit. For finite k and large A^ this is 
clearly useful because the W^'s are defined independently of A^ and p. 

It is not difficult in principle to enumerate normalized /c-plets in standard 
lexicographic order, hence in particular in order of increasing \V\, and then 
compute for each normalized /c-plet the value of \E\. 

We know that 

TFW = N 

and 



Tr Afi = 0. 
For k = 2, the only sequence is (1, 2), so 

WJP = pN^. 

For A; = 3, the only sequence is (1, 2, 3), so 

WJP = p^NK 

For k = 4, the sequences are (1, 2, 1, 2), (1, 2, 1, 3), (1, 2, 3, 2) and (1, 2, 3, 4), 
so 



Tr M4 = p'^m + 2p'm + pm. 

For A; = 5, the sequences are (1, 2, 1, 2, 3), (1, 2, 1, 3, 2), (1, 2, 1, 3, 4), (1, 
2, 3, 1, 2), (1, 2, 3, 1, 3), (1, 2, 3, 1, 4), (1, 2, 3, 2, 3), (1, 2, 3, 2, 4), (1, 2, 3, 
4, 2), (1, 2, 3, 4, 3), and (1, 2, 3, 4, 5), so 



Tr M5 = p^N^ + bp^N^ + bp^N^. 
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For k = 6, one finds 41 sequences leading to 



Tr M6 = p^N^ + {3p^ + Qp^)N^ + {9p^ + Qp^ + 5p^)N^ + {4p^ + Qp^)N^ + pN^ 
For /c > 7, by counting the sequences with \V\ > k — 2, 

TrM^ = /A^^ + {(t - 2A;) / + kp'''^} iV^^ 

+ { (^ _ 5|i + A;2 + 5A;) p'^ + (f - A;2 - 6A;) /-^ + (f + A;) p'^-^} A^^ 

+0(iV*=-3) 

= p^^A^'^ + (-2A;/ + A;/-i)iV'=-i 



Tk-3\ 



With ten days of computation on a workstation, we have obtained the number 
of normahzed A;-plets with given \V\ and \E\ up to /c = 18. The results are 
available upon request to the authors. 

3.2 Sum rule for p = 1 

We have checked our enumeration against a simple sum rule. We put p = 1 
and sum over \E\. In this case, with probability 1, the random graph becomes 
complete and the matrix M is equal to J — Id, where Id is the N xN identity 
matrix and J is the N x N matrix with all entries equal to 1. But J — Id 
has only two eigenvalues, A^ — 1 with multiplicity 1 and —1 with multiplicity 
A^ — 1. So, for p = I, 



Tr M^ = Tr (J - Id)^ = {N - 1)^ + (iV - 1)(-1)^ 
The general formula reduces to 

(AT _ 1)^= + {N- 1)(-1)'= = J2N\n = J2 N^T^Kn 

W n 

where T>k^n is the number of normalized /c-plets W with \V\ = n. Going to 
generating functions, the left-hand side gives 

^N fk 

Y: {{N - ir + {N- l)(-l)^) ^'-r = e-\e-^' + e-(x - 1)), 
while the right-hand side gives 

ryN fk O-fc 



Random incidence matrices: moments of the spectral density 10 

Hence 

k,n '^■ 

By using any symbolic computation software, the computation of all the 
T>k,nS up to, say, fc = 50 takes only a few seconds. We can express Pfc,n 
in terms of standard Stirling numbers of the second kind, Sk^n- We shall 
meet them again in Sec. |5.5|. They c an be characterized by the relation x^ = 
En<Sk,nX-. Using the trick Tr ikf^+i + TWW^ = N{N - 1)'= = J2n<Sk,nN^ 
one finds X'jt+i^n + 1^k,n = <Sk,n~i, which gives for A; > 1 

fc-i 

T^k,n = 2^(~1) ^Sr,n-1- 

r=\ 

For x = 1, we get the generating function of P^ = I]„I^fc,n, the total 
number of normalized fc-plets. That is e^*^ -i-*) ^nd we recognize that the T>k 
are so-called generalized Bell numbers |^loane|| . To give an idea of the size 



of the computations of the first moments, the values of P^ for fc = 0, ■ ■ ■ , 18 
are 1, 0, 1, 1, 4, 11, 41, 162, 715, 3 425, 17 722, 98 253, 580 317, 3 633 280, 
24 Oil 157, 166 888 165, 1 216 070 380, 9 264 071 767 and 73 600 798 037. With 
our combinatorial interpretation of P^, it is clear that P^ ^ ^!- On the other 
hand, the saddle point evaluation shows that logP^ c< k{\o^k + o(logA;)), 
confirming that the growth of the computation is extremely rapid. 

3.3 Model with random signs 

The above calculations could be adapted to the variant of the model (see 
Sec. |1]) defined by a parameter a G [0, 1] where the non-zero elements of the 
random matrix are +1 with probability a and —1 with probability 1 — a. 
Equivalently the edges of the random graph are dressed with a random sign 
±1. 

For an admissible A;-plet / = (zi, ^2, ■ ■ ■ , "^fc), we find Mi^i^Mi^i^ ■ ■ ■ Mi^,i^ = 
p^^{pby° = p%^° where 6 = 2a — 1 is the sign asymmetry, / = /g + Zq is the 
number of distinct pairs among {11,12}, {^2, "^3}, ■ ■ ■ , {ik, h} and le (resp. Iq) 
the number of those which are repeated an even (resp. odd) number of times. 
So, Eq. (0) becomes 

TrM'= = ^A^^p'6'°, (3) 

w 

and the first moments can be exactly computed by enumerating all the nor- 
malized fc-plets: 

WW = pN\ 
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WW = p%^N^, 
WW = p%'^N^ + 2p^N^ + pN\ 
WW = p%^N^ + 5p%^N^ + 5p%^N^, 
WW = p%^N^+{3p%^ + 6p^b'^)N^ + 

{9p^b^ + 6p%'^ + 5p^)N^ + (V + Qp'^)N^ + pN-. 

Of course, the case a = 1 gives previous results. If we concentrate on the 
even model with a = 1/2 (for which 6 = 0), we see that the summation over 
random signs keeps only the walks for which all the edges are visited an even 
number of times. Consequently for all the odd moments, Tr M^'^+^ = and 
the density of states is a symmetric distribution. 

3.4 Generating function 

While very convenient for explicit enumeration, the formula in Eq. (^) is not 
always convenient for theoretical arguments. So we reformulate it. 

Starting from a normalized fc-plet W, we have defined two sets V and 
E. Recall that V is of the form {!,■■■ ,n} for some n and that E is made 
of pairs of distinct elements of V. These are exactly the data for a labeled 
graph with vertex set V and edge set E. In this framework, W can be 
interpreted as a closed walk on the graph visiting all edges (this implies in 
particular that the graph is connected), the order of first visit to a vertex 
respecting the natural order : W starts at vertex 1, and its first visit to 
vertex 2 occurs before its first visit to vertex 3 and so on. Note that in this 
formulation, many labeled graphs do not appear (for example, the labeled 
graphs for which vertices 1 and 2 are not linked by an edge). To resume, only 
some labeled graphs and some walks visiting all edges appear. On the other 
hand, take an unlabeled connected graph G isomorphic to the one defined 
by the normalized sequence W. Clearly, it is possible to label it in such a 
way that W describes a closed walk on G visiting all edges. This labeling 
can be achieved exactly in |ylMt(G)| (the order of the automorphism group 
of G, see Sec. ||) ways. Indeed, two distinct labelings have to describe a 
non-trivial automorphism of G because W determines completely a labeled 
graph isomorphic to G. Written symbolically, this means that 

w G \^ut{G)\ ^^^^ 

where on the left-hand side the summation is over normalized fc-plets whereas 
on the right-hand side the summation on G is over isomorphism classes of 
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connected graphs and the summation on W{G) is over walks on G of length 
k visiting all edges of G. So, the fundamental identity can be rewritten as 

where the sum over G is over isomorphism classes of connected graphs, 
Aut{G) is the automorphism group of G, V{G) and E{G) are respectively 
the vertex set and the edge set of G and Wk{G) is the number of closed walks 
of k steps on G visiting all edges of G. Note that graphs with |V^(G')| > A^ or 
|ii^(G)| > k do not contribute. Although Eq. (^ was established for A; > 2, it 
is valid for k > 0. For A; = 1, as a walk with only one step cannot be closed, 
Wi{G) = and Tr M = 0. For /c = 0, as a walk with zero steps is closed and 
visits one vertex and zero edges, Wo{G) = for every graph except for •, the 
graph with one vertex. For this graph, Wq{») = 1, leading to Tr M^ = N. 
By convention, the empty graph is not counted as connected, so it does not 
appear in Eq. (^). 

The above formula can be used to build a generating function by summing 
over k. On the left-hand side J2k>o -^^Tr Af^ = Tr (1 — XM)~^ is a rational 
function of A because it is the average of a finite number (namely 2^*^^^^'/^) 
of rational functions. To deal with the right-hand side, define Wg{X) = 
J2k ^''Wk{G). This counts closed walks on G of arbitrary length visiting all 
edges of G. Were it not for the last constraint, life would be easy because the 
generating function for the closed walks on G of arbitrary length is simply 
Tr (1 — XG)~^ (in this formula and in later algebraic expressions involving 
graphs and traces, we use a convenient abuse of notation : G denotes at the 
same time the graph and the incidence matrix obtained by labeling it, any 
choice of labeling leads to the same traces). To suppress walks that do not 
visit all edges, we can use inclusion-exclusion to obtain: 

"'-W-^Tz^-g^TrlcaT-H'gTv^-i^.... (5) 

where X^gw is the sum over all subgraphs of G obtained by deleting / edges, 
so this formula ends at / = \E{G)\. This expresses Wg{X) as a finite sum of 
rational functions, so Wg{X) is a rational function. To summarize, we have 
proved an identity between rational functions (so that it is possible to assign 
values to A) : 



Tr i = V i A^n^M pl^(G)liyg(A). (6) 

1-AM ^\Aut{G)\ ^ ^^ ^ ^' 
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4 Fixed edge probability p 



In the large A^ limit with p fixed, Fiiredi and Komlos ||Fiiredi81| (following 



work by Wigner |[Wigner5^ , Arnold [|Arnold67|| and Juhasz |[Juhasz8T|| ) have 



given a detailed description of the spectrum : it consists A^ — 1 "small" eigen- 



values which after a rescaling by a factor Jp{l ~ p)N build up a semicircle 
distribution of radius 2 and one large eigenvalue (the Perron eigenvalue) 
whose distribution is Gaussian with average pN + (1 — 2p) and finite vari- 
ance 2p(l — p). The appearance of this isolated eigenvalue is due to the 
non-vanishing average of the matrix elements. 

We have used this theorem to check our numerical simulations and our 
algorithm for computing the moments. We have also computed numerically 
level spacings. They obey the Gaussian Orthogonal Ensemble statistics with 
good accuracy. 

4.1 Monte-Carlo simulations 

Monte-Carlo simulations consist in generating a lot of incidence matrices of 
random graphs, computing and studying their spectra. 

A matrix is obtained with the following procedure. It is a symmetric 
matrix M of size N x N. Its diagonal elements are set to 0. Its non-diagonal 
elements Mij = Mji (with i ^ j) represent the edges of the graph. They are 
randomly and independently chosen: their values are 1 with probability p 
and with probability 1 — p, where p is a fixed parameter. 

As the matrix is symmetric, it is diagonalizable and all its eigenvalues 
are real numbers. They are computed using the appropriate routine in Nag 
library. This procedure is repeated with several random matrices. As we are 
interested by the asymptotic behavior when A^ is large, we compare different 
sizes of matrix: N= 1000, 2000 and 4000, with 40, 20 and 10 matrices 
respectively. So, for each case, the statistical study is done over 40000 random 
eigenvalues. 

When the edge probability p is fixed, in the spectrum of a given matrix, 
we must make the distinction between its largest eigenvalue and the A^ — 1 
others. Indeed, all the elements are non- negative {Mij > 0). Moreover, 
when A^ is not too small, the random graph is connected and the matrix 
is irreducible. Then, the Perron-Frobenius theorem [pantmacher64|| assures 



that the eigenvalue with the largest modulus (the Perron eigenvalue) is non- 
degenerate, positive, and that the elements of the corresponding eigenvector 
are all positive. 

Numerical observations show with great accuracy that the average of 
this Perron eigenvalue is pN, plus a finite correction which depends on p. 
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Figure 1: Histograms of spectra of random graph incidence matrices, for 
several sizes N of matrices and several edge parameters p. The x-axis is 



rescaled by Jp{l — p)N and the area of each histogram is normalized. All 
the curves coincide with the semicircle of radius 2, drawn with a solid curve. 

Furthermore its variance (i.e the square of the width of its distribution) is 
also finite et depends on p. Indeed, for large N, the Perron eigenvector is 
equal to (1 . . . 1)"^, plus small fluctuations, and the eigenvalue is about pN. 
As the rest of the spectrum has a large A^ behavior which is different we 
eliminate, in the rest of this section, the largest eigenvalue of each matrix. 
On Fig. in histograms of eigenvalues for several values of p are shown. To 



allow comparison, the eigenvalues (x-axis) have been divided by \/p{l — p)N, 
and the y-coordinates have been scaled in order to normalize the area of 
each histogram. We see clearly that the asymptotic shape of the rescaled 
distribution is a semicircle of radius 2. In particular, only the variance (or 
the width) of the distribution depends on A^ and p, but the shape remains 
the same. So already for A^ ~ 1000 the agreement with the large A^ estimates 
in the Fiiredi-Komlos theorem is very good. 

The nearest neighbor spacing distribution is commonly studied to observe 
fluctuations in random spectra. If (Ai, A2, ■ ■ ■ , Aat) are the eigenvalues in 
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Figure 2: Histograms of normalized nearest neighbor spacings in spectra of 
random graph incidence matrices, for several sizes N of matrices and several 
edge parameter p. All the curves coincide with the GOE spacing distribution 
distorted by the semicircle law, drawn with a solid curve. 



ascending order of a given random incidence matrix, we define the normalized 
spacings as 

in order to have (s) = 1, by omitting the last spacing sj^^i which involves the 
Perron eigenvalue. On Fig. 0, histograms of normalized spacings are shown. 
We see clearly that they have the same shape. This indicates that the large 
N behavior is independent of p. 

To allow the comparison with the spacing distribution of the Gaussian 
Orthogonal Ensemble (GOE) of random matrices, we use the "Wigner sur- 
mise" 

9o(s) = -sexp(^--s 

It is considered to be an excellent approximation of the GOE distribu- 
tion ||Mehta9l| ]. However, a direct comparison between the Monte-Carlo 
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histograms and qo{s) would give bad results. Indeed, qo{s) describes the 
distribution of the locally normahzed spacings A^p(Aj)(Ai+i — Aj), where p(A) 
is the eigenvalue probability distribution. So the probability distribution of 

Si = N{Xi+i - Xi) is 

q{s)= f dX p{Xy qo{p{X)s). 



By taking the semicircle distribution p(A) = A/tt^/I — 4A^ with a diameter 1 
in order to have (s) = 1, 

qis) = ^sexpf^sA ^(1,3,1.^) 



12 /-4 2\ ^ r(A: + l/2) 1 /4 

^ . exp(^^s j ^^^ r(l/2) fc!(fc + 2)!U' 



where F(l/2, 3, z) is a generalized hypergeometric function. On Fig. |^, we 
see that the Monte-Carlo simulations coincide with the function q{s), giving 
good evidence that the incidence matrices of random graphs are, as expected, 
in the universality class of GOE. 

4.2 Perturbative expansion for the Perron eigenvalue 

We retrieve quickly the main features of the distribution of the Perron (i.e. 
largest) eigenvalue via a perturbative expansion. By definition, M = p{J — 
Id) where Id is the N x N identity matrix and J is N x N matrix with 
all entries equal to 1, that is, A^ times the projector on \Q) = -4^ (1 ■ ■ ■ 1) 



T 



TV 



So 1^7) is the Perron eigenvector of M with eigenvalue p{N — 1). We define 
D = M — M. Simple manipulations ||Wigner35| show that an eigenvalue A of 



M whose eigenspace is one-dimensional and not orthogonal^ to \Q) satisfies 

l=pN{Q\-—^—-\Q). 
X + p — D 

The perturbative expansion in D gives for the Perron eigenvalue 

x = p{N-i) + {n\D\n) + -\- ({n\D^\n) - {n\D\nf) + ■ ■ 

pi\ ^ ' 

Explicit computation yields 



{Q\D\^) = 0, 



,2 



{n\D^\n) - {n\D\ny = p{i-p){n-i){n-2)/n. 



^ This explain why perturbation theory in D singles out a particular eigenvalue, the 
Perron eigenvalue, the only eigenvalue not orthogonal to jfi) to order in D. 
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Hence 

A^ = 

This is enough to show that the distribution of the Perron eigenvalue has 
average A = pN + 1 — 2p + 0{1/N) and finite width A^ — A = 0(1). Let us 
call tr' M'^ the fcth moment of the distribution of other eigenvalues, 




tr' M'' = — ( Tr iVf^ - A' 



Comparison of our formulae for Tr M*^ and A'^ leads to 



tFM = -p + 2p/N + 0{l/N'^), 

tZW = p{l-p)N + p{3p-2) + 0{l/N), 

tFW = -3/(1 -p)iV + 0(1), 

t7M4 = 2[p{l-p)fN'^ + 0{N), 

ti'M^ = 0{N''^^) for A; > 5. 



This is of course consistent with [[FiirediS 1|| . Note however that for the Lapla- 



cian matrix (see Appendix ^ the hypotheses of the theorem are not fulfilled 
because the diagonal elements are correlated to the rest of the matrix and 
have a variance of order A^. And indeed, the Laplacian has an entirely dif- 
ferent spectral distribution. 

4.3 Comparison with the model with random signs 

The random sign model defined in Sec. |1| (where Mij = +1 with probability 
ap, —1 with probability (1 — a)p and with probability 1—p) is also covered 
by the results in [[burediSll . If a 7^ 1/2, the "small" eigenvalues build a 
semicircle of radius 2 after rescahng by \/p{l — h'^p)N (where, as before, 
h = 2a — 1), and the large eigenvalue has a Gaussian distribution with average 
(A^ — l)bp + (1 — b'^p)/b and variance 2p{l — b'^p). 

As another check of our formulae, we give a short proof of the semicircle 
distribution for the symmetric (a = 1/2) random sign model, where Mij = ±1 
with probability p/2 and with probability 1 — p. To compute Tr M'^ from 
Eq. (^), we must keep in the summation only the W^s for which the edges are 
repeated an even number of times. Then for the odd moments, Tr M2'^+^ = 0. 
For even moments, Tr M^'^, the maximal number of distinct edges in W is 
k. Moreover in the large A^ limit with p fixed, we retain the W^s with the 
maximal number of vertices. This maximum is k + 1 and is obtained with 
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ly's associated to rooted planar trees with k edges, as explained later in 



Sec. |5.5| . Then, for a = 1/2 

1 



Tr M2fc = Ckp'^N'' + 0{N''-^) 

where Ck are Catalan numbers, Ck = {2k)\/[k\{k + 1)!]. As the Catalan 
numbers are the moments of the semicircle law of radius 2, this shows that 
the density of states is the semicircle law of radius 2^JpN . 

5 Fixed average connectivity a 

In this section, we present results for the large A^ limit with a = pN fixed. 
After displaying numerical observations, we prove that the density of states 
has an infinity of delta peaks for any a. Then we explain how to compute the 
2/c-th moment of the density of states which is a polynomial in a of degree 
k. Finally, we give bounds for the coefficients of these polynomials and use 
the bounds to show that the spectrum is determined by the moments. 

5.1 Monte-Carlo simulations and observations 

In Sec. m, we have defined a variant of the model where the non-zero elements 
of the random matrix have a random sign: +1 with probability a and —1 
with probability 1 — a. We will prove in Sec. |5.3| that in the large A^ limit with 
a fixed, the moments of the density of states are independent of a. Then it is 
expected that the density of states evaluated by our Monte-Carlo simulations 
are very similar for those equivalently obtained ||Evangelou92| , [Evangelou92a 



for the even model with a = 1/2. However as we will clarify some points, we 
have repeated these simulations. 

Monte-Carlo simulations have been done with the same procedure as in 



Sec. ^A\ . We have seen that, for a fixed edge probability p, the shape of the 
distribution of eigenvalues is a semicircle. But when p becomes of order 1/A^, 
the distribution gets strongly distorted. To compare different sizes A^, let us 
fix a = pN. It is the average number of I's in a given row (or column) of 
the random matrix. For the graph, it is the averaged connectivity (i.e. the 
average number of neighbors of a given vertex). 

On Fig. ^ histograms of eigenvalues for several values of a are shown. 
For each value of a, three sizes of matrices, A^ = 1000, 2000 and 4000, have 
been simulated with 40, 20 and 10 matrices respectively. So, for each case, 
the statistical study is done over 40000 random eigenvalues. For fixed a, the 
three curves are superposed, in the limits of Monte-Carlo fiuctuations. So 
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Figure 3: Histograms of spectra of random graph incidence matrices, for 
different sizes N of matrices and probability parameter p = a/N, with a fixed. 
The X-axis is rescaled by ^Ja and the area of each histogram is normalized. 
For comparison, the semicircle is drawn for a = 8. 
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Figure 4: Cumulative histograms of random graph incidence matrices, for 
different sizes N of matrices and probability parameter p = a/N, with a = 1. 
Each vertical step represents a delta peak. As the three curves are quite 
similar, they are shifted along the x-axis to be more visible. The largest step 
is at X = 0. 



we can consider that we observe the asymptotic distribution dpa for large A^, 
which depends only on a. To allow comparison between different a's, the 
eigenvalues (x-axis) have been divided by ^/a, and the y-coordinates have 
been scaled in order to normalize the area of each histogram. 

When a is small, we see a forest of delta peaks. This has been previ- 
ously observed in other models of sparse random matrices [[Kirkpatrick72|, 



[Evangelou83|| . Their heights are not representative with this kind of his- 
togram because they depend on the width of the bin, chosen arbitrarily. To 
give a correct representation of the importance of delta peaks, the cumula- 
tive distribution function (i.e. the integral of the distribution from — oo) is 
better. We plot it for a = 1 on Fig ^: each vertical step corresponds to a 
delta peak, with a weight equal to the height of the step. As the heights are 
comparable between the different sizes, the delta peaks survive in the limit 
N -^ oo. 
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We observe that for any a the bigger delta peaks are, in order of impor- 



tance, at X = 0, ±1, ±V2, ±(v^ ± l)/2 (golden mean), ±y/3, ±^2 ± \/2, 
etc. These values can be recognized as eigenvalues of small trees. 

When a increases, the heights of delta peaks decreases, but their positions 
along the x-axis do not move (before the ^/a rescaling). With this kind 
of histogram, when a delta peak becomes too small, it seems to disappear 
because it is drowned in the rest of the distribution. But we will show that the 
spectrum has an infinity of delta peaks for any a. For large values of a, the 
shape is close to the semicircle, obtained when p is finite (which corresponds 
to the limit a = oo). 

In Sec. 0, we have seen that the topology of random graphs has a per- 
colation transition at a = 1. But, this transition seems without effects for 
the distribution of eigenvalues: on Fig. |], distribution for a = 0.8, 1 and 1.2 
are qualitatively similar. In particular, we have checked that the height of 
the delta peak at x = is regular around a = 1. The same remarks apply 
to the localization transition, conjectured [pi)vangelou92| , |Evangelou92a|| to be 
for a ^ lA. 

In contrast, we observe a change of behavior between a = 2.6 and 2.8, 
for the distribution dpa in the vicinity of x = 0. For a < 2.6 (resp. > 2.8), 
dpa{X) decreases (resp. increases) when A goes to 0"^. It is difficult to say if 
the limit dpa{X) when A goes to 0"*" is or not in the small a phase, and cxo or 
not in the large a phase. This transition seems to be related to a transition 
at a = e, where the second derivative of the height of delta peak at x = as 



function of a is discontinuous ||BauerOOa[] . Unfortunately, this discontinuity 



is too small to be seen in our Monte-Carlo simulations. 

We have also studied the distribution of nearest neighbors spacings. Our 
conclusions are similar to those of Evangelou and Economou |[Evangelou92|, 



[Evangelou92a |: for small a, the distribution numerically coincides with an 



exponential, and for large a with the GOE law. In the vicinity of the lo- 
calization transition, a ~ 1.4, the distribution interpolates between these 
two forms. On the other hand, we have not been able to reproduce re- 
sults ||Evangelou92|, [Evangelou92aj| concerning the singularity of the spectral 



distribution as A goes to 0. 

5.2 Existence of delta peaks 

In this section, we explain that, for any value of a, the distribution has an 
infinite number of delta peaks. More precisely, we show that they are delta 
peaks at all eigenvalues of finite trees. However their heights are exponen- 
tially decreasing functions of a; most of them are hidden in simulations by 
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the statistical noise and the distribution seems to be quite smooth for large 
a. 

The delta peaks at tree eigenvalues have, at least, two origins ||Kirkpatrick72 



the small connected components, which are trees, and small trees grafted on 
the giant component (percolation cluster). 

The random incidence matrix is block-diagonal, with one block per each 
connected component. As shown by Eq. |l], the average number of times 
a given tree T appears as a connected component of the random graph is 
proportional to A^. So, for any eigenvalue of the incidence matrix of any 
tree with n vertices, a contribution of height a"'~^e~'"°' / \Aut{T)\ to the cor- 
responding delta peak appears. These eigenvalues are algebric numbers, i.e. 
solutions of a polynomial equation with integer coefficients of degree at most 
n. The height (but not the position) of the peak depends on a. The height 
decreases exponentially with n, so only eigenvalues of small trees appear with 
repetitions in Monte-Carlo simulations. 

Furthermore, for a > 1, even the giant component contributes to delta 
peaks. To explain why, we first define grafting. Let G = (V, E) be a graph, 
G' = (V, E') and G" = {V" , E") be two subgraphs of G, and V^ be a subset 
of V. We say that G is obtained by grafting G' on G" along Vq if Vq = V'nV" 
and E is the disjoint union of E' and E" . 

If this is the case, suppose moreover that the incidence matrix of G' has 
an eigenvector cj)' (with eigenvalue A) whose components on Vq are 0. Then 
the vector obtained by extending cj)' to l^ by is an eigenvector of G with 
the same eigenvalue A. 

Now if G is the giant component of a random graph, if G' has a finite 
number of vertices and is connected, then G' is a tree with probability 1. 



Indeed, it is known |Pollobas98|| that finite connected subgraphs with loops 
are suppressed by powers of A^~^, even if they belong to the giant component. 
For a given tree T, the average number of times G can be obtained by grafting 
a subgraph G' isomorphic to T on a subgraph G" along n points is of order 
Na\^^'^^^e^'^~'^^'^^'^°' /\Aut{T)\ for large a. So we have shown that if T has an 
eigenvector (with eigenvalue A) vanishing on some vertices, A appears as a 
delta peak in the spectrum of the giant component. 

Now if A is an eigenvalue of the incidence matrix of a tree T' with eigenvec- 
tor 0', then A also appears as an eigenvalue of some tree T whose associated 
eigenvector vanishes on some vertices of T: for instance, choose a vertex v on 
T', take two copies T| and T'_ of T' and built a tree T by joining v+ and V- 
to a new vertex vq. Then = (0', 0, -0') on ^(r) = V{T'^) U {vq} U V{T'_) 
is an eigenvector of T with eigenvalue A and this one can be grafted along 
Vq (see Fig. |]). So we have shown that the giant component contributes to 
delta peaks in the spectrum at all eigenvalues of finite trees. 
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Figure 5: A symmetric graft on the giant component gives delta peaks in the 
distribution of eigenvalues. 

To summarize, at any eigenvalue of any finite tree, a delta peak appears in 
the distribution, with contributions both from small connected components 
and from the giant component (for a > 1). We do not know if delta peaks 
appear at other positions, due to other mechanisms. On the other hand, 
we have not been able to prove that for a > 1 the giant component gives a 
continuous part in the distribution of eigenvalues. 

5.3 General considerations 

From an analytical viewpoint, we use equation (|^) which we rewrite as 

w 

and recall that on the right-hand side, E and V can be interpreted as vertices 
and edges of a connected graph. It is then an elementary topological fact 
that \E\ > \V\ + 1, with equality if and only if the graph is a tree. For fixed 
k, the possible graphs form a finite set, and there is no difficulty to take the 
large N limit : 

ak= lim — Tr Mfc = yal^l 

N^^ N V 

where the sum of over normalized fc-plets W associated to trees. The exis- 
tence of the N ^ oo limit above is a strong indication of the existence of 
a limit eigenvalue probability density dpa, for which /ifc is just the usual fc*'' 
moment. 

In the same way, in the large A^ limit with fixed a, Eq. (^ becomes 



1 „ 1 ^ 1 



lim -Tr — - = y -----al^(^)llVy(A), (7) 

N^ocN 1-AM "^ \Aut{T)\ ^ ^' ^^ 

where the sum runs over the trees T. For trees, we observe that the computa- 
tion of Wt{\) using Eq. (^) can be simplified: the sum over all subgraphs of 
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Figure 6: Decomposition of T 



T can be reduced [[BauerOO|] to the sum over all subtrees obtained by deleting 
leaves of T. 

The dominance of trees here and in Sec. |l] has a similar origin. Though 
the giant component contains of order A^ loops, only a finite number of them 
are finite : as already noticed before, a finite connected induced subgraph of 
a random graph at fixed a is a tree with probability 1. 

On a tree, it takes an even number of steps to make a closed walk, and 
/ifc = for odd k. This elementary observation implies in fact that trees 
can be bicolored, and by standard argument, this shows that trees have a 
symmetric spectrum (see Sec. E^ for an application of this property). As 



random graphs look locally like trees, it is not too surprising that they also 
have a symmetric spectrum (this is of course true only for A^ — > oo). Later 
we shall show that the moments fik determine the distribution. Then /i^ = 
for odd k implies that the distribution of eigenvalues is indeed symmetric. 
For the even moments, fi2k is a polynomial of a with degree k, 

I 

where Ik^i is the number of normalized 2fc-plets W associated to trees with / 
edges. 

If we consider the variant of the model with random signs (see Sec. |l|), the 
same arguments apply. In particular, only VT's associated to trees contribute 
to the sums. But in this case, any edge is visited an even number of times 
and, following Eq. (^, /i2fc does not depend to a. Then, in the large A^ limit 
with a fixed, the parameter a is irrelevant. 



5.4 Recursion relation for Xi 



kj 



Consider a normalized 2k-]ilet W associated to a tree T, with k > 1. Recall 
that W induces a walk covering T, hence a labeling of T (vertices are labeled 
in order of their appearance in the walk). Let r be the root of T, the vertex 
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labeled 1, where the walk starts. There is always an edge between vertices 
labeled 1 and 2. If the edge {1,2} is cut, the tree breaks in two trees, T' 
whose root r' is vertex 2 and T" whose root r" is vertex 1 (note that r" = r). 
The tree T can be seen as trees T' and T" linked by the edge {r",r'}, as 
shown in Fig. |^. Note that T' and T" are arbitrary trees. They could for 
instance consist of a single vertex. 

The walk W can be decomposed accordingly : the walker starts at vertex 
r", walks along edge {r", r'}, makes a closed walk on T' (made possibly of zero 
steps), walks along edge {r',r"}, makes a closed walk on T" (made possibly 
of zero steps), and so on and so forth, and finally comes back to vertex r", 
after covering T' and T". Let 2k' (resp. 2k") be the number of steps made 
on T' (resp. T"). We can glue the small pieces of walks on T' (resp. T") 
together into a single walk, covering T' (resp. T") because W covers T. This 
walk is a sequence of vertices on T' (resp. T"), and the "label and substitute 
algorithm" applied to this sequence leads to a normalized 2A;'-plet W (resp. 
to a normalized 2/c"-plet W"). 

Now, the 2/c-plets W giving rise to the same pair {W, W") are easily 
counted. They differ from each other by the organization of the steps on 
edge {r',r"}. There are 2m = 2{k — k' — k") such steps. If m! (resp. m") 
is the number of returns of W (resp. W") at vertex r' (resp. r"), there are 

exactly 1 (resp. 1 ) possibilities to insert the m 

steps from r' to r" (resp. the m steps from r" to r'). Note that the —1 in the 
above counting comes from the fact that the last (resp. first) step from r' to 
r" (resp. from r" to r') is fixed. The total number of visits at vertex r" is 
m + ?7i", the first term counting visits from the edge {r' ,r"} and the second 
visits from T" . If V (resp. T") has /' (resp. /") edges, T has / = /' + /" + 1 
edges. To summarize, the number X^ ^ „ of normalized 2fc-plets associated to 
trees with / edges and containing m times the number 1 satisfies the following 
recursion relation, for A; > 1 

T _Y^-r T fm' + m 

-Lk,l,m — Z^-Lk',l',m' -Lk",l",m" 




where the sums runs over non-negative indices A;', k'\ I', I", in\ m" and rn 
with relations k' + k" + m = k, I' + I" = I — 1 and m + m" = ni. 

Note that Tk^i^m vanishes ii I > k (every edge is visited at least twice 
for a covering closed walk on a tree), if m > A; (two successive terms in a 
normalized 2/c-plet are distinct, so 1 cannot appear more than k times) and 
if 771 = unless A; = / = (every non void normalized 2A;-plet contains 1). 
Finally, Tqqq = 1. This gives more than enough boundary conditions to 
compute the X^ ^ ^ recursively. 
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3248 


3731 


2002 


429 








8 




254 


3628 


14062 


23020 


18928 


8008 


1430 






9 




510 


10805 


57516 


127029 


144024 


91392 


31824 


4862 




10 




1022 


31740 


227030 


654395 


968544 


828495 


426360 


125970 


16796 



Table 1: The number of normalized 2k-plets associated to trees with I edges. 

On a workstation, a day of symbolic computation with this formula gives 
the Xfc,/,m as integers up to A; = 50. With 12 digits precision, a Fortran 
program goes to k = 120 in about the same time. The results are available 
upon request to the authors. 

Note that the index m is not directly relevant for the computation of 
moments, because X^.z = J2m^k,i,m, but we have not been able to obtain a 
closed recursion without the index m. Coefficients Xk,i for small k and / are 
given in Table |1|. 

Let us note that each moment is a polynomial in a, so that in particular it 
does not exhibit any singularity at the percolation or localization transition. 
This is not very surprising, because by the general algorithm (any A^ and p) 
/ifc is computed by exploring connected subgraphs with at most k sites of the 
random graph. But as noted before a finite connected induced subgraph of 
a random graph at fixed a is a tree with probability 1. To see any trace of 
a transition, one would have to look at the behavior of fi2k for large k. For 
instance, we tried to see if complex zeroes of ^2k in the variable a (or other 
related quantities) have a tendency to accumulate near the real axis. We 
have found no conclusive evidence. The rapid growth of the coefficients Tk,i 
makes a numerical study up to A; = 50 difficult, even if we know all numbers 
exactly. 



5.5 Bounds for Xi 



ki 



In this section, we will show that X^^u the number of normalized 2fc-plets 
associated to trees with / edges satisfies the bounds 

Sk,i < 1k,i < Ci Sk^i, (9) 
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where Ci are Catalan numbers defined by 

(2/)! 



C, 



11(1 + 1] 



and Sk^i are Stirling numbers of the second kind (i.e. the number of ways 
of partitioning a set of k elements into / non-empty subsets) which can be 
computed with the formula 

5m = I E (-1)'-"' (;„) ."'. (10) 

m=0 \ / 

Let W he a. normalized 2k-plet associated to a tree T. We view W as a. 
walk on T. Then W allows to put more structure on T. First, the starting 
point of the walk turns T into a rooted tree. This allows to talk about sons 
of a vertex v, the root being the initial ancestor. Its sons are its neighbors, 
and so on. Then W also endows the sons of a vertex with an ordering : the 
order of first visit. This means that W naturally endows T with a structure 
of plane rooted tree, with the convention that the root is at the top of the 
tree, and the elder son is always the leftmost one. Then, 

Ik,i = / .IkjTi) 

Ti 

where the sum runs over the plane rooted trees Ti with / edges, and 1k{Ti) is 
the number of admissible walks on T; with 2k steps, where admissible means 
starting and finishing at the root of T; and visiting all the I + 1 vertices by 
respecting the order of birth among brothers (i.e. a vertex can be visited 
only if its brothers on its left have been visited before). 

First, we will prove that Ik{T*) = S^^i for a particular tree — the star-like 
tree — which gives the lower bound of Eq. @. Then, we will prove that 

Ik{Ti) <Sk,i (11) 

for any tree T/. As the number of plane rooted trees with / edges is Ci (see 
e.g. [^tanley97|| ), this gives the upper bound of Eq. (^. 



The star-like tree T* is made of a root r and / sons: {si, S2, . . . , s^}, 
see Fig. ^. If / = 0, then I^^Tq) = 6kfi = S^fl. So we consider now 
that / > 1. All the normalized 2k-plet associated to T* can be written 
as (r, s„(i), r, s„(2), ■ ■ ■ ,r, Sjj{k)) where v is an onto map from [1, k] to [1, /], be- 
cause an admissible walk is made of k successive double-steps from r to a son 
and return, and visits all the sons. The inverse map v~^ makes a partition of 
[1, k] into / non-empty subsets: v~^{j) = {i\l < i < k; v{i) = j} for j G [1, /]. 
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Figure 7: Labels of vertices of the "star" tree T*. 

Conversely, such a partition gives a admissible walk, because the targets j 
associated to each subset are uniquely determined by the birth rule, with the 
following process: the subset containing 1 is labeled by j = 1 and removed, 
the subset containing the smallest remaining number is labeled by j = 2 and 
removed, etc, up to the last subset labeled by j = /. This is a one-to-one 
correspondence between the admissible walks and the partitions of [1, k] into 
/ non-empty subsets, which are counted by Stirling numbers of second kind. 
Hence, 2fc(T*) = Sk,i for the star-like tree T^. 

Let us now consider a given plane rooted tree Ti with / edges. So Eq. ([TT|) 
is true for k < I because 1k{Ti) = (the walk is too short to visit all the 
edges) and ior k = I because 1i{Ti) = 1 = Si^i (an admissible walk on T/ of 
length 21 exists and is completely determined by the birth rule). Note that 
this implies that T;/ = Ci. For k > I, we will prove Eq. ( [Tl| ) by induction, 
assuming Eq. (|ri|) to be true for every {k',l') when k' < k. 

Let r be the root of T; and r' the leftmost son of r. As shown on Fig. |^, 
we break Ti in three parts, the edge {r, r'}, the sub-tree T' with root r' and /' 
edges and the rest T", which is a sub-tree with root r" = r and /" edges. We 
have I = I' + I" + 1. An admissible walk W on Ti with 2k steps is composed 
of 2k' steps on T', 2k" steps on T" and 2m steps on the edge {r,r'}, with 
k = k' + k" + m and m > 1 . 

If T is a plane rooted tree, we call 1k,m.{T) the number of admissible walks 
on T with 2k steps and m returns to its root, and define 

m=l \ / 

The arguments used to establish Eq. (H) can be repeated to show that 

rh+k'+k"=k 
1k{Ti) = 22 'Hk' ,rh{T')Hk" ,m{T") . (12) 

m,k'k" 



Now 



k' _ 

'Hk'^miT') < ^ lk',m'iT')\ - _. =1k'iT' 



m'=l 
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ai a2 a/. 



r =r 




bj b2 b/" 



Figure 8: Labels of vertices of the "bi-star" tree TiV 



I" 



and by the induction hypothesis, the last term is at most 

k' + m — 1\ 



Ik'iTi,] 



m 



liTii) because for each walk on T^t, 



ni 



k'. Hence, 



But this is precisely Tifc'.m 

'^k',miT') < 'Hk',miTii). As the same argument holds for T", 2k{Ti) < 

Ik{T^},) where T^},, is the "bi-star" tree with T' = T*, and T" = Tf„. 

It remains to show that Xfc(T;**„) < Sk^v+vi+i- If /' = 0, it is true because 
the bi-star is simply the star T*. As Eq. (|12D is symmetric by exchange 
between T' and T", it is also true for /" = because Xfc(^z*'o) ~ ^fc(^o'l') — 
Sk,i- 

So we have just to consider the case I' > 1 and /" > 1. The vertices of 
the bi-star are labeled as shown on Fig. p. The admissible walks are divided 
in two classes: (I) the walks finishing by a step on the right subtree T*,, and 
(II) the walks finishing by a step from r' to r" . 

For a walk W of class (I), we call W'^ the walk W without the two last 
steps; W^ is made of 2fc — 2 steps. The class (I) is divided in two sub-classes: 
(Ii) W^ has not visited all the vertices. Then W = (W\r",bi") and W^ is 
an admissible walk on T^*„_i. (I2) W"^ has visited all the vertices. Then 
W"^ is an admissible walk on T^*n and there is /" choices to built W from 
lyl". In total, the number of walks of class (I) is Ik-i(Tiiy,^^) + /"Jfc_i(T;**„), 
bounded above by Sk-id-i + l"Sk-i,i. 

For a walk W of class (II), we call W'^ the walk W with the first and last 
steps removed. Then W'^ is made of 2A; — 2 steps, starting and finishing at r' 
(and not r) and covering the bi-star. The delicate point is the moment of its 
first visit to r. We have /' + 1 sub-classes. In sub-class (Ho), the first visit of 
r is before the first visit of ai. In sub-class (IIj) with 1 < i < /' — 1, the first 
visit of r is after the first visit of Oj and before the the first visit of flj+i. In 
sub-class (II/O) the first visit of r is after the first visit of a;/. By using the 
birth rule, each sub-class corresponds to the admissible walks on one of the 
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M I'-i 



i = 





i = l' 



Figure 9: The V + 1 plane rooted trees corresponding to the class (II). 

plane rooted trees with / edges drawn on Fig. |^. For each of these /' + 1 trees, 
the number of steps is 2k— 2, so that the induction hypothesis apphes. Hence, 
for class (II), the total number of walks is bounded above by (/' + l)Sk^i,i- 
Then, Ik{T^},) < Sk-i,i-i + (/" + /' + 1)5^-1,/ = Sk-i,i-i + lSk-i,i. But, the 
Stirling numbers obey to the recursion relation Sk,i = Sk-i,i-i + lSk-i,i, so 
that we have reached our goal: Eq. (|TTD is proved and Eq. (^) follows. 

We can now establish two important features of the eigenvalue distribu- 
tion dpa '■ it's support is unbounded, but it Fourier-Laplace transform is an 
entire function, and in particular, the distribution of eigenvalues is deter- 
mined by its moments. 

By using the property of Stirling numbers n^ = I]JLo <Sk,in- and by sum- 
ming it with X^^^^/n!, one obtains 



Sk{x)=J2Sk,ix' = e-^Yl 



oo ^k 



-X" 



1=0 



n=0 



n\ 



This leads to a crude bound for the Stirling polynomials when x is real 
positive. In fact, Sq{x) = 1, Si{x) = x, S2{x) = x^ + x and for A; > 3, 
Sk{x) < k^ + ga^'(fc-i) because for k >?> and any n, n^ < k^ + fc". 

To see the unboundedness of the support of the eigenvalue distribution, we 
observe that the Stirling polynomials iSfc(a) are the moments of even order of 
the even a-dependent probability measure ^ S^o ^('^(2^~v^)+'^(^+"\/^)) 
on the real line. The support of this measure is clearly unbounded. As the 
moments of dpa are larger than the Stirling polynomials, the support of dpa 
has to be unbounded. 

To see the properties of the Fourier-Laplace transform, we just show 
that term by term expansion of J^^ dpa{X)e'^^ in powers of s leads to a 
series with infinite radius of convergence. For this, we just need to bound 



/-r^P.(A) 



(2fc)! 



P2k/{'^k)\. We know that this is at most CkSk{a)/{2k)\ 



fe!(fc+l)!'-''=v'^) 



Ski a). We have seen than for A; > 3 this is less than 



fcfc_)_ga(fc-l) 

A:!{A;+1)! 



which 
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for large k behaves like , k(k+iv ~ ^^^ ^'^ ^°^ ^^^ ^' ^°' ^^ principle, 
the knowledge of the moments determines the probability distribution dpa- 
However, it is not easy to extract accurate local information on dpa from the 
knowledge of a finite number of moments. 

5.6 Special values 

Eq. (H) can be used to compute X^^i for special cases. This is best done using 
generating functions. Instead of giving the details, let us just note that this 
leads to closed forms for X^ / for fixed / and any k or for fixed k — I and any /, 
leading to a a satisfactory description of the borders of the table of moments. 
For k = 0, Iqi = 6qi. For A; > 1, the first cases are : 



Xfe,o 


= 0, 


Ik,i 


= 1, 


Ik,2 


= 2'^-2, 


Xfc,3 


= 3''-^ + Lj'' + LU'' 



3 ■ 2^' + 2 



with uj + uj = 3 and ujuj = 1. 
On the other border 

III = -TTTZ TTT = Cl. 



^1+1,1 

^1+2,1 
^1+3,1 



11(1 + 1) 

(2/ + 2)! 
(/-!)!(/ + 3)!' 

(2/ + 4)! /2^11/ + 2 
(/-!)!(/ + 6)! 2 ' 

(2/ + 6)! /^ + 32/3 + 323/2 + 232/ - 48 
(/-!)!(/ + 9)! 6 ■ 



That Ii^i is just the Catalan numbers is not surprising for two reasons. 
First combinatorially, for any plane rooted tree with / edges, Ti{Ti) = 1, so 
Ii,i = Cl, the number of plane rooted trees. Second, we know that for fixed 
p, the distribution of eigenvalues is governed by the semicircle law. It is 
not surprising that when a goes to infinity, the same distribution reappears. 
Indeed, the fact that Tii = Ci is equivalent to the fact that the typical size of 
eigenvalues is y/a, and that after the rescaling A = y/ax, dpa{\) converges 
to semicircle law ^\/4 — x'^dx^ for which the even moments are the Catalan 
numbers. 
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The above equations are reminiscent of the meanders problem ||Dil:'rancesco96 | 
In particular, we observe that diagonals of Table |l] verify 

i2l + 2u)\ P2u-2il) 

'+"'' (/-!)!(/ + 3m)! u\ 

where P2m-2(0 is a polynomial with integer coefficients and with leading term 
^2«-2_ gy conjecturing this form, the coefficients could be determined for the 
first values of u, from the exact knowledge of the first Ii-^u,i- Unfortunately, 
we know no general formula for I^^i. 

The asymptotics of Xk,i for fixed / > 2 and large k are governed by a 
simple relation. 



As explained in Sec. ^75t Ik,i = J2Ti^kiTi), and Xfc(T/) is maximal for the 
star tree T;* as for the bi-star tree T**'^^ q (isomorphic to the star but with a 
leaf as root), with Xk(T*) = Xk(T**^Q) = Sk,i. Hence, Xk,i > 2Sk,i for / > 2. 
But, for / fixed, we will show that I^^i is asymptotic to this lower bound when 
k is large, because among all the Ci plane rooted trees, the contributions of 
T* and T**^ q become dominant. Eq. ( p!0D says that Sk,i ~ l''/U for large k. 

Let M be the incidence matrix of a given plane rooted tree T; with / edges: 
the number of admissible walks with 2k steps, Ik{Ti), is bounded above by 
Tr M^'^ which counts the closed walks on T with 2k steps. So Ti contributes 
for large k only if the Perron- Frobenius eigenvalue of M is not smaller than 

Vi. 

A tree can be bicolored: if we use black (B) and white (W) as the colors, 
each vertex is either black or white, in such a way that edges connect only 
vertices with opposite colors. Then, for an eigenvector with eigenvalue 
A, the vector 0', defined by 0^ = 0, on black vertices and — 0j on white 
vertices, is eigenvector associated to eigenvalue —A. Hence, the spectrum is 
symmetric. 

As Ti has / edges, TrM^ = J2i,j M^j = 21, and for the spectrum of M, 
X]j A^ = 21. We see that T; contribute only if its spectrum is { — \//, 0,v^} 
where is / — 2 times degenerated. By noting $ the Perron- Frobenius eigen- 
vector, associated to y/l, and $' the eigenvector associated to — V^, in this 
case 

As $' is obtained by inverting the signs of $ on white vertices, Mij = when 
i and j have the same color and Mij = 2vZ$j$j when i and j have different 
colors. The Perron theorem assures that $j > for all i. Hence, Mij = 1 
when i and j have different colors. As loops are forbidden in a tree, one of 
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the colors must color only one vertex. So M must be the incidence matrix 
of a tree isomorphic to the star T^, which prove Eq. (pTSl). 

Using the saddle point approximation, it is not difficult to show that for 
large k and / with k/l fixed, Sk,i grows faster than any exponential of /. In 
this regime, logXfc,/ ~ logiS^,; because logC/ ~ /log4 <^ logSk^i- But the 
asymptotics of 2k,i seem much more complicated. 

6 Conclusion 

In this paper we have made a detailed study of the spectral density of large 
random graphs incidence matrices, both in the fixed edge probability p and 
in the fixed average connectivity a limits. 
For fixed p, our results are: 

• A simple general algorithm to compute arbitrary moments (polyno- 
mials in N and p), easy to implement on a computer (but of almost 
factorial growth), leading to an explicit form for the ffist 18 moments. 

• Numerical diagonalizations of Monte Carlo samples of random matrices 
and analytical arguments to show that when the random signs are not 
chosen with a symmetric probability, the main change is that a "Perron 
Frobenius" eigenvalue has to be eliminated to recover the semicircle law 
and the GOE results. 

For fixed a (the finite connectivity limit), our main contributions are: 

• Numerical diagonalizations of Monte Carlo samples of random matrices 
for different values of a, with a curious observation for a ~ 2.7. 

• A general proof that the spectrum contains delta peaks at all eigen- 
values of tree incidence matrices, and that, slightly surprisingly, they 
receive non- vanishing contributions even from the infinite cluster when 
a > 1. 

• A recursion relation of combinatorial origin to compute the moments. 
Odd moments vanish, and we have explicitly computed the first 120 
even moments with the help of a computer. Our algorithm counts 
objects which are related to many other tree enumeration problems of 
independent interest. 

• A proof that the growth of the moment is slow enough so that they 
determine the spectral distribution entirely. 
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• The moments are insensitive to the random signs, and so is the spec- 
trum by the previous remark. 

• The moments are polynomials in a. In particular, they do not exhibit 
any singularity at the classical or quantum percolation transitions. 

• However, we have given a general formula (|^) to prepare the ground 
for a more refined study pauerOOa|| of the delta peak at the eigenvalue 



A = in the spectral distribution. This delta peak is directly relevant 
to quantum percolation and is non-analytic in a at a = e, which we 
believe is connected to the curious numerical observation at a ~ 2.7 
mentioned above. 

One of the natural continuations of this work would be a careful numerical 
and/or analytical study of moments of large order, to see if the percolation 
transitions have a measurable impact on global characteristics of the spectral 
density. 

From a more mathematical point of view, and because of their close con- 
nection with delta peaks in the spectra of random graphs, it might be of 
interest to be able to characterize the eigenvalues of tree incidence matrices 
of a given size and how their distribution and spacings evolve for large sizes. 

Acknowledgments We thank one of the referees for pointing reference 
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A Random Laplacian matrix 

In this appendix, we give results for the Laplacian matrix on random graphs. 
These are obtained using an adaptation of methods previously described. 

The random Laplacian matrix L of size N is defined as L = Z^ — M 
where M is a random incidence matrix (see Sec. |^) with vanishing diagonal 
elements and D is the diagonal matrix whose element Da = J2j=/=i ^ij is the 
connectivity of the vertex i in the random graph associated to M. So all 
rows or column sums of L vanish. 

The Perron-Frobenius theorem applies to N — L: the Perron eigenvector 
is the uniform eigenvector (1, ■ ■ ■ , 1)"^, associated to the eigenvalue Ap = 
for L. It follows that the spectrum of L is real and non-negative. Moreover 
the multiplicity of is the number of connected components. 



To compute Tr L'', the method explained in Sec. ^]l] is adapted to the 
Laplacian. After expanding (D — M)'^, and eventually using the invariance of 
the trace by cyclic permutation, we must compute several terms individually. 
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Of course, the diagonal elements Da need a special treatment. With the 
trick Da = Y^j^i^ij = Hj^i^ij^ji^ Da becomes a double-step {i,j,i) in 
the description as a walk on a graph. Now, an admissible fc-plet is made up 
of single-steps (for M) and double-steps (for D). After enumeration, 

WL = pN^, 

TFT2 = p^N^ + 2pN\ 

THJ = p^N^+{<op^ -p^)N^ + ApN\ 

WT? = p'^N^+{12p^-3p'^)N^+{25p'^-Qp^)N^ + 8pN^. 

A computer program can expand these results for a dozen of k. They can be 
checked with sum rules: for p = 1 the matrix is deterministic and Tr L'' = 
{N - 1)N''. 

For finite A^ there are in general several connected components, hence the 
Perron eigenvalue Ap = appears with multiplicity. However for fixed p the 
random graph is connected in the large A^ limit, and Ap = appears only 
once. So it is reasonable to consider centered moments nik = (A — A)'^ where 
the means run on the other A^ — 1 eigenvalues. From previous equations, 

r7i2 = 2p(l — p)N, 

ms = Ap{l-p){l-2p)N, 

m4 = p{l-p)[p{l-p){9N -A2) + 8]N. 

In the large A^ limit with p fixed, ma/rTig goes to and m^/ni^ goes 
to 9/4. The limit distribution has a bell-like shape [|Biroli99| , Pavagna"99 



intermediate between the semicircle and the Gauss law for which iriii/m^ is 
2 and 3 respectively. 

In the large A^ limit with a = pN fixed, the /c-plets contributing to the 



dominant term of order A^ are associated to trees, as described in Sec. |0 . 
So the kth moment is a polynomial in a of degree k, 



Ilk = lim — Tr /^^ = XI ^k,ia\ 

where £fc,/ is the number of normalized fc-plets (with single or double steps) 
associated to trees with / edges. Note that now k is not constrained to be 
even, due to the double steps. Following the arguments and notations of 
Sec. |5.4|, we call Ck,i,m the number of normalized /c-plets associated to trees 



with / edges and containing m, times the number 1 (i.e. with m, return to the 
root of the associated tree). 

In comparison with Sec. |5^ , as the steps on the first edge {r',r"} are 
single or double, we call u' (resp. u") the number of such steps finishing on r' 
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k\l 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


1 


1 




















2 


2 


1 


















3 


4 


6 


1 
















4 


8 


25 


12 


1 














5 


16 


90 


85 


20 


1 












6 


32 


301 


476 


215 


30 


1 










7 


64 


966 


2345 


1722 


455 


42 


1 








8 


128 


3025 


10696 


11659 


4928 


854 


56 


1 






9 


256 


9330 


46453 


71082 


43779 


12012 


1470 


72 


1 




10 


512 


28501 


195340 


404540 


342642 


135357 


26040 


2370 


90 


1 



Table 2: The number of normalized 2k-plets, with single or double steps, 
associated to trees with I edges. 

(resp. r"), i.e. single steps {r",r') and double steps {r',r",r') (resp. {r',r") 
and {r",r',r")). Then for k>l, Ck,i,m satisfies the recursion relation 



^k,l,m — / ,^k\l',m'-^k'\l" 



m' + u' 



U 



m" + u'' 
u" - 




where the sum runs over non-negative indices fc', k'\ I', I", m', m", u' and u" 
with relations k' + k" + u' + u" = k, I' + 1" = I — 1 and u" + m" = m, with the 
convention that the first binomial coefficient must be taken as 1 (and not 0) 
when m' = u' = 0. To start the recursion relation, we need the boundary 
conditions for A; = 0: Co,i,m = ^i^m- 

By summing on m, we retrieve the coefficients Ck,i = J2m ^k,i,m- The 
first of these are given in Table 0. As the average /ii = a is not zero, we 
compute the centered moments m^, which are still polynomial in a of degree 
the integer part of k/2. The coefficients Ci^^i are given in Table 0. The shape 
of the distribution evolves with a: in particular, it becomes non symmetric as 
the odd centered moments do not vanish. Again, we observe that the results 
for large N with fixed p are found by keeping the last diagonal of Table § 
and Table ^, corresponding to the large a limit. 
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